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We use the concept of a KirchhofT resistor network (alternatively random walk on a net- 
work) to probe connected graphs and produce symmetry revealing canonical labelings of the 
graph(s) nodes and edges. 

1 Voltages 

Electrical networks and random walks on networks have been the object of study 
for a while (see Doyle, Snell [2]). Recently these concepts have been taken up by 
data miners and network analysts to extract useful information from data and 
networks [3] [4] [6] . Before this the same techniques where used by computational 
chemist to label molecules [5] . If one takes Kirchhoff 's Current Law as a given 
the exposition of what takes place on a network is easier to follow than using 
random walks. Given an undirected connected graph G — (V,E), with vertex 
set V = {1, 2, ...N} and edge set E = {1, 2, ...M} the graph can be represented 
by an NxN adjacency matrix A where Oy = w ii there is an edge of weight 
u; > between nodes i and j otherwise = 0. Make the edge weights 1 ohm 
resistors. Choose one node a to be a current source / and one node 6 to be a 
current sink — /. Now we apply KCL. 

iy = for i ^ a,b j connected to i 
Y.j iaj = I 

iij — Cij{vi — Vj) Cij = conductance 
16{k - a) - 15{k -b) = J^j Ckj{vk - Vj) 

For one ohm resistors Cy = aij aij e A where A is the adjacency matrix. 
Let D represent the diagonal matrix of vertex degrees Da = . 
Let Ci = (0, 1, ...,0n) 
Let 1=1 amp then 



(Ca - 66) = {D- A)V 



= Lv 

Where L is the graph's Laplacian. L is singular and not invertable. There 
is one undetermined degree of freedom. We can solve the KCL equations by 
using the Laplacian's pseudoinverse Z/+ or we can remove one vertex from the 
Laplacian and invert the resulting N — 1 matrix setting the removed vertex 
voltage to zero. 

V = L+{ea - Cb) + Xe 

Where Xe is an undetermined constant. 

We may also add a universal sink node (an edge to every other node) to the 
graph to get a A''a;A'' ivertable matrix. Here we can disregard the sink node. 

2 Comparisons 

We now determine the node voltages v over all possible node pairs. This yields 
a set of A'' voltage vectors with Nx{N — 1) elements were Nx{N — l)/2 of 
the elements are the negative value of the others. These voltage vectors may 
be subtracted from each other to yield a set of M current vectors. To see if 
two graphs are isomorphic one does a comparison of the M current vectors of 
the graphs. The vector elements may be rearranged at this point, to aid the 
comparison as their order is immaterial. To reveal a graph's symetries compare 
the A'' node voltage vectors to each other (element order is immaterial) and place 
them into groups of identical vectors. These groups of nodes are the orbits of 
the automorphism group of the graph. 
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